ACEEE Int. J. on Communication, Vol. 02, No. 03 , Nov 20 1 1 



QRC-ESPR1T Method for Wideband Signals 



Nizar Tayem 

Electrical Engineering Department 

Prince Mohammad Bin Fahd University 

Al Khobar 3 1952, Saudi Arabia 

ntayem@pmu.edu. sa 



Abstract - In this paper, a new algorithm for a high resolution 
Direction Of Arrival (DOA) estimation method for multiple 
wideband signals is proposed. The proposed method proceeds 
in two steps. In the first step, the received signals data is 
decomposed in a Toeplitz form using the first-order statistics. 
In the second step, The QR decomposition is applied on the 
constructed Toeplitz matrix. Compared with existing schemes, 
the proposed scheme provides several advantages. First, it 
requires computing the triangular matrix R or the orthogonal 
matrix Q to find the DOA; these matrices can be computed 
with 0(n 2 ) operation. However, most of the existing schemes 
required eignvalue decomposition (EVD) for the covariance 
matrix or singular value decomposition (SVD) for the data 
matrix; using EVD or SVD requires much more complex 
computational 0(n 3 ) operation. Second, the proposed scheme 
is more suitable for high-speed communication since it 
requires first-order statistics and a single snapshot. Third, 
the proposed scheme can estimate the correlated wideband 
signals without using spatial smoothing techniques; whereas, 
already-existing schemes do not. Accuracy of the proposed 
wideband DOA estimation method is evaluated through 
computer simulation in comparison with a conventional 
method. 

Index terms - array signal processing, uniform linear array, 
Wideband signal, Toeplitz matrix. 

I. INTRODUCTION 

Source localization using a sensor array is one of the key 
techniques in mobile communication. It plays an important 
role in many application fields such as radar, sonar, mobile 
communication, Multiple-Input Multiple-Output (MIMO) 
system, biomedical, and wireless sensor networks [1], [2]. 
Numerous classical subspace methods e.g., Multiple Signal 
Classification (MUSIC) [3], Estimation of Signal Parameters 
via Rotational Invariance Techniques (ESPRIT) [4], and 
Maximum Likelihoods method (ML) [5] have been developed 
to estimate the DOA of the narrowband signals. Many of 
these narrowband methods are not applicable to wideband 
signal directly. The reason is that the energies of narrowband 
signals are concentrated in the frequency band that is 
relatively small compared with center frequency. The 
wideband signals encountered in many applications such as 
high-speed data communication, passive sonars, and speech 
signal processing . A wide range of techniques have been 
proposed to deal with the wideband source localization 
problem. The maximum likelihood extended for the wideband 
signal scenario in [6] and [7] provides optimal statistic 
properties, but it requires a multi-dimensional search and is 
highly non-linear. The one dimensional search method in [8] 
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proposed for non-coherent signals requires less 
computational load. A class of representative high -resolution 
DOA estimation approaches for wideband sources when the 
signals are highly correlated due to multipath or reflection is 
the coherent signal subspace method (CSSM) in [9] and [10]. 
The CSSM converts the wideband signal subspaces into a 
predefined narrowband subspace by focusing matrices; 
subsequently, the narrowband subspace-based DOA 
estimation methods, e.g. MUSIC, can be applied. However, 
the CSSM method requires initial DO As, and its performance 
is sensitive to these initial values. Two other classes of 
wideband source localization approaches are based on 
orthogonality of subspace. In one class the test of 
orthogonality of projected subspace (TOPS) [1 1] is referred 
to as a non-coherent method. The method requires singular 
value decomposition (SVD) to get the signal subspace of 
each frequency point. However, using the SVD required the 
computational cost of 0(n 3 ) operations. In the other class, 
the test of orthogonality of frequency subspace (TOFS) [12] 
has been proposed for non-coherent wideband signals . TOFS 
constructs the searching steering vector of every possible 
DOA and frequency. However, TOFS cannot resolve the 
desired DO As when the signal to noise ratio (SNR) is low. 
The objective of this paper is to provide an efficient algorithm 
to estimate DOAs for multiple wideband sources and to 
reduce the computational complexity when compared with 
the existing schemes. To achieve these objectives, the 
proposed method proceeds in two steps. In the first step, the 
received signals data is decomposed in a Toeplitz form using 
the first-order statistics. In the second step, the QR (no 
previous reference for an acronym) decomposition is applied 
on the constructed Toeplitz matrix. Compared with existing 
schemes, the proposed scheme provides several advantages. 
First, it requires computing the triangular matrix R or the 
orthogonal matrix Q to find the DOA. The computation of 
QRfor ftXft with the Toeplitz matrix [13-14] requires 0(n 2 ); 
whereas, the existing schemes require eigenvalues 
decomposition (EVD) or singular value decomposition (SVD) 
for an n x n cross-spectral matrix. The major drawback of 
SVD or EVD is that they require much more computational 
complexity 0(n 3 ). In addition, the proposed method is more 
suitable for high-speed communication first order statistics, 
a single snapshot, and requires computing the triangular 
matrix R or Q only, which further reduces the computational 
complexity. However, the existing method involves multi- 
dimensional research or high order statistics which increases 
the complexity and computational load even more. In addition, 
the proposed method is able to resolve the DOAs of highly 
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correlated sources. This paper is organized as follows: in 
section II, the proposed method is introduced based on first 
and second order statistics. In section III, we present the 
simulation results. A conclusion is given in Section IV 

E. SYSTEM MODEL AND PROPOSED METHOD 

Consider a symmetric uniform linear array (ULA) 
consisting of M=2N+1 identical omni-directional antennas, 
where the distance between adjacent antennas is d. Consider 

L wideband signal sources from the directions l ,0 2 ,---,0i. 
Then the received signal by the m-th sensor can be expressed 

as 



X,n(t)=Y J S,{t-T m {e i )) + n m {t) 



(1), 



where sf^t) is the Z-th source signal; n (i) is the additive white 
Gaussian Noise (AWGN) at the m-th sensor; and 

r m [Of ) = (m)d cos 0, 1 v is the delay, where v and d denote the 
propagation speed of source signals and the distance 
between sensors, respectively. A narrowband DOA 
estimation and source localization problem requires that the 
inter-sensor spacing should satisfy the so-called half- 
wavelength constraint to avoid the phase ambiguity caused 
by the multi-valued property of the cosine function. However, 
half-wavelength constraint is not valid to a multiple wideband 
source data model. Applying the Discrete Fourier Transform 
(DFT) to (1), we obtain 



i-i 



-njat) m=^ r m,-M,-£ (2). 



If we use the element at the center of the array as a reference 
point, then the [2N + l)x 1 output vector can be written as 

x(a> k )=[x_ N (co k ),---, x (co k ),---, x N (co k )J (3), 

where the superscript (Y represents transpose. We can 
express X(co k ) as 



x (<°k)= I>/(<»,tM®i >#/)+« W 



/=i 



(4), 



where a{w,6[) is the array response vector in spatial- 
frequency domain and is defined as 

a{(0 k ,8) = [e ]COtT -«^ ••• 1 ••• e 1C0tT ^ e) ] (5). 

The matrix formulation of (4) can be written as 

X(a k )= A(o) k )s((D k )+n((o k ) (6), 

where ^)=W{^\'-^iM7 an M^)=["-^V"%k)F 
are the DTFT data vectors of the source and noise, 

respectively. The [2N + l)x L matrix 



A(o) k ) = [a(co k ,0 l ),---,a(co k ,0 L )] 



(7) 



is the array manifold matrix in the spatial-frequency domain. 
Equation (6) means that in the time-frequency domain, the 



model of wideband and nonstationary source localization is 
similar to the one of narrowband and stationary case at each 
time and frequency bin. For the first step in the proposed 

method, we convert the output data vector x(co k ) into a 
Toeplitz Hermitian data matrix with 
dimensions [N xl)x[N xl). The advantage of introducing 
the Toeplitz Hermitian data matrix is that it has a rank that is 
related to the DOA of the sources irrespective of whether the 
sources are coherent or not. This matrix, denoted Z, is given 
by 



X-nM 

x -(N-l) Wt I 
x -(N-2) V^k I 



>¥ 



x (co k ) 


x -iM 


x-iM 


x iM 


x oM 


X-i(co k ) 


x 2 M 


x iM 


x (co k ) 



x N (co k ) x N _ x (co k ) x N _ 2 (co k ) 



iM 



(8) 



The signal space of the Toeplitz data matrix in (8) has full rank 
which implies that all the incident sources can be detected 
whether they are coherent or non-coherent. The second step 
of the proposed method is to calculate the R data matrix or Q 
data matrix from the received Toeplitz data matrix (8) to estimate 
the DOA for form coherent multiple wideband sources. By 
using the Toeplitz data matrix in (8) the necessary information 
about the noise subspace or the signal subspace can be 
extracted using QR factorization. One of the reasons that QR 
factorization [15-17] is widely used in adaptive applications 
is that in RRQR the signal information can be effectively 
updated making the algorithm suitable for tracking moving 
sources. In addition, the computation of QR for the nxn 
Toeplitz matrix [13-14] requires 0(n 2 ) operations. In the 
proposed method, we calculate the triangular factorization R 
to find the DOA for the wideband signal. Using the shift 
invariance property of Toeplitz matrix we can partition the 
data matrix *}> in two ways: 



¥ 



or 



T 



x (co k ) 


f 
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F-i. 




F-i 


K 
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KT 

z x 


)K) 





(9) 



(10), 



where F_ { is an (n-l)x(n-l) submatrix of 4> , 
/ = k-i(fijt) ■ ■ ■ JtM«c)] = an > ,iT = [x_ v {e> k ) 






Let R be the upper triangular factor of the Cholesky 
factorization of 4> r 4> , 



R T R = w Ty ¥ 



(11), 
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where R is the upper triangular matrix from the QR 
decomposition of 4> . The matrix R can be partitioned in two 
different ways 



R 



'1,1 >fr 
O R h 



(12) 



or 



'\c 



r N,N 

U ' 



(13), 
O r N 

wh er e r fi = [r 12 ■■■ r UN J , r lc =[r lN ■■■ r N _ 1-N J ; R b is 

the bottom submatrix and R t is the top submatrix have a 
dimension (N — 1) x (N -1) . From equations (9), (11) and (12) 



we get 



Similarly, from equations (10),(ll)and(13)we get 



1,1 r \ffr 

7/l,l r fr r fr + K^b 



rJr, Rjr lc 

Tn T.J- 

r k / V 1c r \c +r n,n 



T T 1 T_ ] +z K z KT 
y RT T -1+4^) yy 



z K *o{(Ok)+T r y R 

RT M4%J 



(14). 



(15) 



Comparing the left upper submatrices in (15) we get 

RjR, = T T l T_ l +z K z Kr (16); 

while comparison of the lower right matrices on both sides of 
(14) gives us 

r] r r fr + RlR h = yy T + t\t_ x (17). 

Finally, (16) and (17) give the main relation 

Rj,R b = RjR, + yy T - z K z KT - rj r r fr (18) 

where from (14), Ay is given by 
x oMy T+ z T ' r -i 



r fr 



w^f 



(19) 



Equations (18) and (19) represent the base in calculating the 
matrix R. For more detailed information on obtaining the 
elementsofRmatrixreturnto[13] and [14]. After the matrix R 
is calculated, the DOA of the wideband signal can be found, 
as discussed in the following paragraph. Now the calculated 

data matrix R with dimensions (jV + l)x(7V + l)can be 
partitioned as 



R 



R n R, 
R 



12 



22. 



(20), 



where R n is (LxL); R [2 is Lx(n + 1- L) ; and R 12 is 
(N - L + l)x (N - L + 1) . The QR factorization in (34) is called 
a rank-revealing QR factorization if R22 has a small norm. 



Since Rj2 is very small, the basis of the noise space can be 

obtained from the above R factor as follows. Let V be a 
permutation matrix which represents the row and column 

interchange. Then, given the small norm of the matrix R22 , 
the upper triangular matrix R can be written as 

R = [R U R I2 ]V T (21), 

where the matrix R in (2 1 ) defines the null space of Y. Let G 
represent any vector in the null space of R, i.e. RG — • To 

find the structure of G, we partition G into g { with L 

components and g 2 with (N-L+l) components. Then GR=0 
implies that 
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[ R u R n] 



-0 



(22) 



so that Ri\g\ +R1282 = 0- Since 7?^ is a non-singular 
matrix, g x can be written in terms of g 2 as follows 



g x =-R n R n g 2 

Then G can be written as 



Si 



■Rn R, 



(23). 



(24). 



To find the basis of null space R, we choose any set of N-L+l 
linearly independent vectors; for example, the columns of 
the N-L+l dimensional identity matrix. The basis for the null 
space of the upper triangular matrix R, which is also the null 
space of Y, is therefore given by 



H, 



-R n R 12 



1 N-L+l 



(25). 



It is important to observe that the columns of the null space 
of H q r are not orthonormal. This is in contrast to the null 

space which can be derived from the S VD or EVD techniques. 
We now note that the subspace spanned by the columns of 

H qr is orthogonal to the subspace spanned by the columns 

of B, where the columns of B contain the information about 
the AOAs of incident sources. This is similar to the well- 
known MUSIC algorithm in which the eigenvector of the noise 
subspace is orthogonal to the steering vector of the signal 
subspace. To find the AOAs, we search for the minimum 

peaks of # ?r B(#) . Since the basis of is not orthonormal, 

we use the orthogonal projection onto this subspace in order 
to improve the performance by making the basis of the null 
space of orthonormal. So the columns of the null space of 
become orthonormal as follows 



H = H qr \H q r H qr r H qr 



(26). 
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We now employ the root MUSIC algorithm to find the 
wideband DOAs of the incident sources. The proposed QRC- 
ESPRIT using root MUSIC [18]can be written as 



P~ x = a {& t , 6 f H H^a {s> k , Q) 



(27). 



The roots of the polynomial are obtained from (27) for all the 
frequency bins k, and then they are averaged to derive the 
final DO A of the incident signals 



A k=k, 



(28), 



where \k L , k H ] denotes the index of frequency ranges which 
satisfy K = k H - k L +1 . 

m. SIMULATION RESULTS 

Performance of the proposed method is compared 
with that of the TLS-ESPRIT and Root MMUSIC with spatial 
smoothing, which are considered one of the best DOA 
estimation algorithms. The specifications of the simulation 
are shown in Table I. 

TABLE I. 
SPECIFICATIONS OF SIMULATION 



Arrav craifi Eviration 


1 1 -element^unLicrni linear array) 


Ni—faer c: coherent incident 

signals 


2 


DOAs sf incident signals 


30 = and 45'' 


J.CCVul^tlCn tvpe 


CDMA 


D=ti length 


1024 


Sisuil baBdindSi 


:z\:hz 


Cmier Frecuencv 


2 GHz 


SXR rm.se 


0-20 dB 


Snnplmg frequency 


20 MHz 


DFT points 


255 


Frequency raise 


1 KHz to ICMHz 


Xumber of trials 


100 



Figures 1 and 2 show the root-mean-square-error (RMSE) in 
degree for Wideband DOA estimation error versus SNR form 
to 20 dB for two coherent sources at 30° and 45°. It is observed 
that the proposed methods result in lower RMSE compared 
to the TLS-ESPRIT and Root MUSIC with spatial smoothing. 
In addition, the proposed method provides low complexity 
and OOmpiiati onal OOSt 0(n 2 ) compared to the method in [1 1] 
which requires EVD or SVD and 0(n 3 ) operation. Furthermore, 
the proposed method does not require spatial smoothing 
techniques to estimate the coherent sources; whereas, TLS- 
ESPRIT and Root MUSIC do. In Table II, the computational 
cost of the proposed method is evaluated in comparison with 
TLS-ESPRIT and Root-MUSIC methods. Note that the unit 
of computation time is in seconds. We observe that the 
proposed method requires smaller computational cost 
compared with TLS-ESPRIT and Root-MUSIC. 

TABLE II. 
COMPUTIONAL COST 



Data 
length 


TLS-ESPRIT with 
spatial smoothing 


Root MUSIC with 
spatial smoothing 


Proposed 
Method 


256 


0030 


0.050 


0015 


512 


0.055 


0.095 


031 


1024 


0104 


0.173 


0.053 



rv. CONCLUSION 

A new approach is introduced for the estimation of DOA 
for wideband signals. The observation output is arranged 
into a Toeplitz matrix which enables us to perform the 
estimation using one or more snapshots and estimates the 
coherent source without any additional processing compared 
to TLS-ESPRIT and Root MUSIC that need spatial smoothing 
techniques. The proposed methods employ QR factorization 
to estimate an R matrix or Q matrix which has a much lower 
complexity and computational cost 0(n 2 ) compared to existing 
methods, which require EVD or SVD and 0(n 3 ) operation. 
Through simulations we have shown that the proposed 
method exhibit significant performance improvement over 
many subarray methods such as TLS-ESPRIT and Root 
MUSIC. These advantages make our proposed method 
appropriate for real-time implementation. 
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Figure 2. RMSE for wideband DOA angle estimation errors versus 
SNR for Source 2 at 45° by using the proposed method, Root 

MUSIC, and TLS-ESPRIT algorithms. 



Figure 1. RMSE for wideband DOA angle estimation errors versus 

SNR for Source 1 at 30° by using the proposed method, Root 

MUSIC, and TLS-ESPRIT algorithms 
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